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^ , Abstract: We consider high energy colhsions of two shock waves in AdSs as a model of ultra- 
<D ', relativistic nucleus-nucleus collisions in the boundary theory. We first calculate the graviton field 



produced in the collisions in the NLO and NNLO approximations, corresponding to three- and 
four-graviton exchanges with the shock waves. We then consider the asymmetric limit where the 
K>" [ energy density in one shock wave is much higher than in the other one. In the boundary theory 
[ this setup corresponds to proton-nucleus collisions, with the nucleus being the denser of the two 
^ ] shock waves and the proton being the less dense one. Employing the eikonal approximation we find 
^ ■ the exact high energy analytic solution for the metric in AdSs for the asymmetric collision of two 
O delta- function shock waves. The solution resums all-order graviton exchanges with the "nucleus" 
^ shock wave and a single-graviton exchange with the "proton" shock wave. Using the holographic 
' renormalization prescription we read off the energy-momentum tensor of the matter produced in 
proton-nucleus collisions. We show in explicit detail that in the boundary theory the proton is 
5-1 ■ completely stopped by strong-coupling interactions with the nucleus, in agreement with our ear- 
■ ■ ■ ' lier results [1]. We also apply the eikonal technique to the asymmetric collision of two unphysical 
delta-prime shock waves, which we introduced in [1] as a means of modeling nuclear collisions with 
weak coupling initial dynamics. We obtain a surprising result that, for delta-prime shock waves, 
the multiple bulk graviton exchange series giving the leading energy-dependent contribution to the 
energy-momentum tensor terminates at the order of two graviton exchanges with the nucleus. 
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1. Introduction 

In this paper we continue our earlier investigation [1] of colliding shock waves in AdSs. Due to 
the Anti-de Sitter space/conformal field theory (AdS/CFT) correspondence [2-5], the problem 
of two colliding gravitational shock waves, while an important problem from the standpoint of 
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gravity theory [6-9] , may also be relevant for high energy hadronic and nuclear collisions at strong 
coupling [10-14]. 

One of the most important problems in the field of ultrarelativistic heavy ion collisions is the one 
of isotropization and thermalization of the produced medium. There is a growing consensus in the 
heavy ion community that the medium produced in heavy ion collisions at RHIC is strongly coupled 
[15-23]. The challenge for the theoretical community is to understand (i) how the medium, which 
is initially very anisotropic with zero or negative longitudinal component of the energy- momentum 
tensor [24-29], evolves into an isotropic medium described by ideal (Bjorken [30]) hydrodynamics, 
and (ii) why this transition happens over extremely short time scale of 0.3 -^ 0.6 fm/c, as required 
by hydrodynamic simulations [15-22]. 

There is also a widespread belief in the community, supported by a broad range of phcnomcno- 
logical evidence, that the very early stages of heavy ion collisions are weakly-coupled, i.e., they are 
described by the physics of Color Glass Condensate (CGC)/parton saturation [24,26,31-47] (for a 
review of CGC see [48-50]). It appears that the system produced in heavy ion collisions evolves 
with time from the weakly-coupled CGC state to the strongly coupled quark-gluon plasma (QGP) 
described by the ideal hydrodynamics. There are two types of transitions that the system has to 
undergo in order for such process to take place. First of all, at some point in time the system should 
undergo a transition from weak coupling to strong coupling. Second of all, at a (presumably) differ- 
ent time the system will evolve from the anisotropic early state, in which transverse and longitudinal 
pressure components in the energy-momentum tensor are drastically different, to the isotropic later 
state, in which all pressure components are equal (or almost equal) to each other, as required by 
the ideal (viscous) hydrodynamics. We will refer to the latter transition as the isotropization tran- 
sition. The isotropization transition is a necessary condition for the thermalization of the produced 
medium. 

In this work we will assume that the isotropization transition takes place after the strong 
coupling transition. Hence we will study the onset of the isotropization in the strongly-coupled 
framework. Since strong coupling dynamics in QCD is prohibitively complicated, especially for the 
ultrarelativistic processes at hand, we will employ AdS/CFT correspondence, assuming that the 
bulk properties of the collisions and the produced medium in A/" = 4 super- Yang-Mills theory are 
not too different from QCD and would allow us to make conclusions which are at least qualitatively 
applicable to the real life. 

Attempts to study isotropization and thermalization in the AdS/CFT framework have been 
made before. A gravity-dual of Bjorken hydrodynamics was constructed in [51-55]. To obtain it the 
authors of [51] assumed that the medium produced in heavy ion collisions is rapidity- independent. 
Imposing a no-singularities requirement [51] (or simply demanding that the metric is real [28]) one 
then obtains the asymptotic late-time geometry corresponding to Bjorken hydrodynamics. However, 
this result by itself does not prove that Bjorken hydrodynamics is a consequence of a heavy ion 
collision. In other words, it is not clear which early-time dynamics (or, in general, which events in 
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the past) lead to this dual-Bjorken geometry. 

To address this problem, by analogy with the perturbative approaches [24,25,37,56], it was 
suggested in [10] that one should study collisions of two shock waves in AdS space: following the 
dynamics of the strongly-coupled medium produced in such collisions one would be able to see 
how the ideal hydrodynamic state is reached by the medium and whether this late-time state is 
rapidity- independent. In [11] the case of shock wave collisions in the 1 + 1 dimensional boundary 
theory was considered and solved exactly in AdSs geometry. Unfortunately the lower dimensionality 
of the problem severely limits the physical behavior of the produced medium, and does not allow 
to formulate the problem of isotropization. The case of realistic 1-1-3 boundary theory was first 
addressed in [12] using AdSs space with the infinitely-thin delta-function shock waves. The authors 
of [12] constructed a perturbative series for the energy-momentum tensor of the produced strongly 
coupled matter. 

In [1] we generalized the results of [12] by solving Einstein equations in a more general frame- 
work, which does not depend on the exact profile of the shock waves, i.e., whether they are delta- 
functions or some other objects with finite extent. We identified the perturbation series of [12] with 
a series in bulk graviton exchanges with two shock waves (see e.g. Fig. ^ below for an example 
of a term contributing to the series). Most importantly, in [1] it was argued that in a collision of 
any two physical shock waves, they stop shortly after the collision, possibly forming a black hole. 
In the boundary theory this behavior corresponds to the colliding nuclei stopping shortly after the 
collision, probably leading to Landau hydrodynamics description of the system [57]. Such complete 
nuclear stopping would lead to complete stopping of the baryon number carried by the nuclei. As 
such a complete baryon stopping is not observed at RHIC (and, in fact, baryon stopping at mid- 
rapidity at RHIC is rather small [58] in accord with perturbative calculations [59,60]), this indicates 
that colliding shock waves may not be adequate for the description of realistic nuclear collisions in 
AdS. Indeed, an AdS description would apply if the collisions were strongly-coupled at all times: as 
the early stages of RHIC heavy ion collisions are weakly coupled, an AdS/CFT description of the 
collision at all times can not be valid. In an attempt to resolve the issue we suggested in [1] that one 
could use unphysical shock waves with the delta-prime profile. Such shock waves appear to have no 
stopping. It is possible that using delta-prime shock waves as external sources for the AdS/CFT 
correspondence would yield a more realistic description of heavy ion collisions, and would allow one 
to tackle the problem of isotropization in the strongly-coupled framework. 

In this paper we further explore shock wave collisions. In Sect. ^ we extend the expansion in 
graviton exchanges from [1] to two higher orders. We calculate the next-to-leading order (NLO) 
and next-to-next-to-leading order (NNLO) corrections to the result of [1] for both delta-function 
and delta-prime shock waves (see Eqs. (|2.30| and (|2.35|) , along with Eq. (|2.39|) ). 

We continue in Sect, ^by constructing the resummation procedure in which graviton exchanges 
with one shock wave are resummed to all orders while the interaction with another shock wave 
is restricted to a single graviton exchange (see Fig. ^ below). The diagrams are analogous to 



- 3 - 



those resummed in the study of classical gluon fields produced in proton-nucleus collisions in the 
perturbative CGC framework [61-65] (see [50] for a review). We apply the eikonal approximation 
to Einstein equations, which allows us to construct an exact solution for the energy-momentum 



(Eikonal approximation in AdS/CFT was studied before in [66-70].) Our solution would receive 
energy-suppressed corrections if shock waves of finite width are considered. We note that the energy- 
momentum tensor (|3.38|) is not that of ideal hydrodynamics, indicating that the system does not 
reach isotropization/thermalization in proton- nucleus approximation to the collision. Resumming 
graviton exchanges with the nucleus shock wave to all orders allows us to demonstrate the stopping 
of the proton shock wave explicitly. The relevant component of the energy-momentum tensor of the 
proton in shown in Eq. (p.45|) : one can explicitly see that it goes to zero as the light cone coordinate 
(in which direction the proton was initially moving) is increasing. 

We also apply the eikonal treatment to the delta-prime shock waves. The results are quite 
interesting: we show that in the eikonal approximation the series in graviton exchanges terminates 
at the level of two graviton exchanges with the nucleus shock wave. Thus the NLO result for the 
energy-momentum tensor is, in fact, exact for the case of proton-nucleus collisions! The energy- 
momentum tensor for delta-prime shock waves is shown in Eq. ( p.71|) . It is clear from Eq. (|3.71|) that 
the produced medium distribution has a strong rapidity dependence. Therefore it seems unlikely 
that rapidity-independent Bjorken hydrodynamics geometry of [51] could result from a collision of 
two shock waves in AdSs space, though indeed a further study of the full nucleus- nucleus scattering 
problem is needed to unambiguously answer this question. 

We will conclude in Sect. ^ by summarizing our main results. 

2. Perturbative Expansion in Graviton Exchanges 
2.1 General Setup 

Consider a collision of two ultrarelativistic nuclei. Assume for simplicity that the nuclei have 
infinite transverse extent and the same longitudinal thickness at all impact parameters. The energy- 
momentum tensors of the two nuclei can be written as (Ti (x^)) and (T2+4.(x^)) with the brackets 

(. . .) denoting the averaging in the nuclear wave functions and the light cone coordinates defined by 
= (x° ± x^)/-\/2 where is the collision axis. The geometry of the collision is shown in Fig. [l|. 
As was argued in [51], the geometry in AdSs dual to each one of the nuclei in the boundary 
theory is given by the following metric 



tensor of the produced medium in the case of delta-function shock waves, given 





2 dx^ dx + ti{x ) z"^ dx 



^ + dx]_ + dz^] 



(2.1) 



for nucleus 1 and by 



ds^ 




2dx^dx + t2{x'^) z'^ dx' 



,+ 



^ + dx\ + dz^] 



(2.2) 



4 



for nucleus 2. Here dx\_ = {dx^y + {dx'^Y with x^ and x^ the transverse dimensions which we 
will denote using Latin indices, e.g. x*. L is the curvature radius of the AdSs space and z is the 
coordinate describing the 5th dimension with the boundary of the AdS space at 2 = 0. We have 
also defined 

2 7r2 



and 



ti{x ) 



111 



(2.3) 
(2.4) 

m 



(2.5) 



in accordance with the prescription of holographic renormalization [71]. The metrics in Eqs. 
and (|2.2| ) are exact solutions of Einstein equations in the empty AdSs space 

4 _ 

Our goal is to construct the geometry in AdSs dual 
to the collision of two shock waves given by Eqs. ( |2.1| ) 
and ( |2.2| ). In [1] we argued that the single shock wave 
metric in Eq. (|2.1|) (or in Eq. ( p.2|) ) corresponds to the 
single-graviton exchange between the source nucleus at 
the boundary and the point in the bulk where the met- 
ric is measured. The solution of Einstein equations 
( p.5|) for the collision of two shock waves can therefore 
be represented as a sum of tree-level graviton exchange 
diagrams, as shown in Fig. 0. There the source nu- 
clei are represented by thick crosses, with nucleus 1 
given by the crosses on the top, and nucleus 2 given by 
the crosses at the bottom. As was argued in [1], each 
rescattering in nucleus 1 brings in a factor of ti{x~) 
into the metric, while each rescattering in nucleus 2 
brings in a factor of t2{x^). The large thin cross in 
Fig. 1^ denotes the point in the bulk in the argument of 
the metric, i.e., the point where the metric is "measured". One encounters similar diagrams but 
with gluons and in 4 dimensions for nuclear collisions in the framework of McLerran-Venugopalan 
(MV) model [32-34], as was worked out in [37,38,41,56,61,72]. 

Inspired by the graviton-exchange analogy of Fig. ^ we write the metric dual to the full collision 
as [1] 




nucleus 1 



nucleus 2 



Figure 1: The space-time picture of the 
ultrarelativistic heavy ion coIUsion in the 
center-of-mass frame. The colhsion axis is la- 
beled x^, the time is x^. 



ds' 



I? 

z 



— { — 2 dx~^ dx + dx\ + dz^ -|- ti I 



X 



z'^ dx + t2{x'^) z'^ dx'^'^ + 0{tit2 



(2.6) 
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Figure 2: Diagrammatic representation of the solution of classical Einstein equations for the collisions of 
two shock waves. The wavy lines denote graviton exchanges between the sources at the boundary (thick 
crosses) and the bulk. The large cross denotes the point in the bulk where one measures the metric. The 
upper row of thick crosses denotes rescatterings in nucleus 1, each of which generates a factor of ti. The 
lower row of thick crosses denotes rescatterings in nucleus 2, each of which generates a factor of t2- 

Indeed the interesting unknown part of the answer is in the term denoted o{ti ^2) in Eq- ( |2.6| ): this 
term comprises all higher order graviton exchanges, i.e., higher powers of ti(x~) and t2{x^). The 
first term in this expansion, the term proportional to ti ^2 was found in [1] . For a particular form of 
ti{x~) and t2{x~^) given by delta-functions, the expansion to several higher orders in ti and ^2 "was 
constructed in [12]. 

To construct a series in graviton exchanges for a general form of ti and ^2 and to set up the 
general problem we write 

ds"^ ~ ~ ] ~ [2 + G{x~^ , x~ , z)] dx~^ dx~ + [ti{x~) z"^ + F{x'^ , x~ , z)'\ dx~'^ 

+ t2(a;+)/ + F(x+,x",2) c?x+^ + [l + i/(x+, x", 2)] c/xi + (2.7) 

The unknown functions F(x+ ),G{x^,x , z) , and H {x~^ , X , 2;) contain all higher 

powers of ti and ^2- Note that as Eq. (pT|) and Eq. (|2.2|) are exact solution of Einstein equations 
(p.5|), the functions F, F, G, and H contain at least one power of ti and t2 each [1]. 

Substituting the metric of Eq. ( p. 71 ) into Einstein equations ( p.5| ) yields a very complicated 
system of non-linear equations. It is likely that the solution of these equations is only possible 
numerically. Here we will build on the results of [1] to construct the first few steps of the perturbative 
expansion: we will construct the next-to-leading order (NLO) and the next-to-next-to-leading order 
(NNLO) corrections to F, F, G and H. NLO corrections resum terms containing tft2 and tit"^, 
while NNLO corrections include terms with tlt2, t\t\ and tit\. (The powers of ti and ^2 should not 
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be taken literally, they only indicate the number of times ti and ^2 enter the expression.) In Section 
^ we will use the eikonal approximation to resum one power of ti and all powers of ^2- 

Before we start the calculations let us first point out that, according to the prescription of 
holographic renormalization [71], if we expand the unknown coefficients of the metric in Eq. (|2T 
into a series in powers of z"^ 



oo oo 



n=0 n=0 

00 00 

G{x+, x-,z) = z^Yl Gn{x^, X-) z^"", H{x+, x-,z) = z^J2 Hn{x+,x-) z^"", (2.8) 

n=0 n=0 

then the expectation value of the energy-momentum tensor of the matter produced in the collision 
in the boundary theory is given by the first coefficients in the expansion in Eq. ( |2.8| ): 

(T^^) = 7r-2 ^o{x\ x~) (T-) = Foix\ X-) 

1 /V^ .. 
= -:.^Go(x+,x-) {T^) = -^^5^^H,{x+,x~). (2.9) 

Einstein equations ( p.5|) impose two constraints on the energy-momentum tensor: tracelessness 

(T/) = (2.10) 

and energy-momentum conservation 

^,{T^''') = 0. (2.11) 

Imposing the constraints ( p^.lOD and ( |2.11| ) on the energy-momentum tensor in Eq. ( |2.9| ) we easily 
see that the energy-momentum tensor can be expressed in terms of a single unknown function: 

(T+-) = -^Hoix-',x-) (r^) = -^5^^Hoix+,x~). (2.12) 

Here we defined the following integrations 

x+ x~ 

^[...](x+)= I rfx'+[...](x'+), ^[...](x-)= I rfx'-[...](x'-). (2.13) 

—00 —00 

Eqs. (|2.12| ) demonstrate that only one metric coefficient in ( |2.7| ) is needed to construct the 
energy-momentum tensor of the produced matter in the boundary theory. 
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2.2 NLO Results 
2.2.1 NLO Calculation 

To systematically include the graviton exchanges of Fig. ^ into the metric of Eq. ( |2.7D we expand 
the coefficients of the metric in powers of ti and t2- We start by writing 



F{x , x~ , z) 
F{x~^, x~ , z) 
G{x'^ , x~ , z) 



X If* 



= H^^\x+,x-,z)+H^^\x-^,x',z) + H^'^\x^,x-,z) + ... 



(2.14) 



where the superscript (0) denotes terms containing ^1^2, i-e., quadratic in t's, the superscript (1) 
denotes terms cubic in t's (i.e., terms containing tlt2 and tit^)? the superscript (2) denotes terms 
quadric in t's, etc. Note that the expansion in t's in Eq. ( p^.llD is independent of the expansion in 
z's in Eq. (|2.8| ): each term in the expansion in Eq. (p.l4|) can in turn be expanded in powers of z"^ 



as was done in Eq. (|2.8| ), and vice versa. 

The leading order (LO) terms in Eq. ( 2.14|) denoted by the superscript ^ 
(0) were found in [1]. For completeness let us quote the results here 

1 



F(°)(x+,x" 




F(°)(x+,a;- 




G^''\x^,x' 




if(°)(x+,x~ 




We defined 


[1] 


hoi^x'^ , x" 


) = 


Xii^x'^ , x' 


) = 



— Ai(x^,x ^z^ 
-A2(x+, X") z"^ 



-dlho{x+,x-)z^ 
6 



— d"^ hiix'^ ,x 



l:dlho{x^,x ) - ^ dlhi{x^ , X )/ 



-2 ho{x^, X') / - 2 hi{x^, x~)z^ + - ti{x~) t2{x^) z 



hoi^x^jX )z'^ + hi{x^,x )z^. 




92 92 
d 



^4 



ti{x )t2{x^), hi{x^,x ) 

9_ 



/io(a; , X ), X2{x~^,x ) 



3 9+9_ 
ho{x^,x') 



(2.15) 

I2 

Figure 3: Graviton 
ti{x )t2{x^) diagram correspond- 
ing to the LO solu- 
(2.16) tion of Einstein equa- 
tions found in [1]. 



The diagram corresponding to the LO solution given by Eqs. ( ^.151 ) and 
( p.l6| ) is shown in Fig. |^. Note that the terms in Eq. (|2.15|) adhere to the 
pattern outhned in Eqs. ( p.l2| ). 

To find the NLO terms denoted by superscript (1) in Eq. (|2.14| ) we substitute the metric (|2^ ) 
with the coefficients expanded according to Eq. ( p.l4|) into Einstein equations (|2.5| ). Expanding the 
Einstein equations to the cubic order in t's yields the following equations for G^^'^ and H^^^ 



+ 67 Z + 6q Z + 611 Z 



11 



zz 



Gi^) + 2iyf ) - z G^^l -2z H^^} + 48 aj z'^ + 80 ag z'^ + 120 an z 



11 



0. 



(2.17a) 
(2.17b) 
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The coefficients ^7, 6g, ^n, ay, ag, an are known functions of ti and ^2 the exact form of which 
is not important here. The subscripts z, and x~ indicate partial derivatives with respect to 
these variables. Eqs. ( p.l7a| ) and ( p.l7b| ) are labeled according to the lowercase Einstein equations 
components. 

Solving Eq. ( |2.17a| ) for G'i^'' and substituting the result into Eq. ( |2.17b| ) yields the following 
equation for H^^^ 

-3 Hi^^ + 3z Hi^] - Hi^], + 2z^ H^^}^_ ^ + 12 + 16 ^9 + 20 t/^h z^^ = (2.18) 

with the coefficients ip given by 

4 

3 

^9 = ^ [t2ix^)ho ^-^-{x^,x') +ti{x')ho ^+^+{x'^,x~)] 
3 

V'li = ^ [t2(a;+)/ii . (2.19) 

To find the solution of Eq. ( p.l8|) we follow the strategy used in [1] . We ffist expand H^^^ into 
a series in powers of z"^ 

00 

H^^\x+,x-,z) = z^Yl H^2\x^.x')z^"- (2.20) 

?i=0 

Substituting Eq. ( |2.2CI| ) into Eq. ( |2.18| ) and requiring that the coefficients at each power of z on the 
left hand side are zero yields the recursion relation 

(-2n) {2 + n) + H^^l^. + ^7 5„,2 + ^9 5n,3 + ^11 5nA = n>l. (2.21) 

Arguing just like in [1] that causality requires the series (|2.2CI| ) to terminate at some finite order, we 
see that the series can only be terminated if H^^ = 0. The solution of Eq. ( |2.18| ) is thus given by 

H'~^\x+,x-,z) = H^^\x+,x-)z^ + H[^\x+,x-)z^ + Hi^\x+,x-)z^ + Hi^\x+,x~)z''^ (2.22) 

with the coefficients 

6 96 , 2880 

^'^--a7^^-(a^^- (2.230 

Hi'^ = --^ijn. (2.23d) 



components of Einstein equations which are not shown here explicitly we can find (and have found) 
F(^) andF(i). The remaining components of Einstein equations do not generate further constraints. 
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Note that G^^\ F^^^ and F^^') are indeed needed to 
construct the metric at higher orders in the expansion 
in t's. However, as we argued above and as shown in 
Eq. (|2.12| ), only Hq^^ is needed to obtain the energy- 
momentum tensor of the produced matter at NLO. 
bmce at NLO G^^\ and F^^^ are not needed for 
the boundary theory physics that we are interested in 
here, we will not present explicit expressions for these 
quantities. 

The NLO solution in Eqs. (|]2|) and (|]23D is 
represented diagrammatically in terms of graviton ex- 
changes in Fig. ^. As shown in Fig. ^ the NLO solution 
consists of a single rescattering in one nucleus and a 

double rescattering in another nucleus. As can be seen from Eqs. ( [^.23D and ( ^.191 ), NLO solution 
includes terms with two powers of ti and one power of ^2 and terms with two powers of t2 and one 
power of tl. 



^2 ^2 ^2 

Figure 4: Graviton diagrams corresponding 
to the NLO solution of Einstein equations 
found in this Section. 



2.2.2 Delta-Function Shock Waves at NLO 

It is instructive to find what the obtained results give for specific shock waves described by particular 
forms of ti{x~) and t2{x^). Define the transverse pressure p of the produced medium by 



Combining Eq. with Eqs. ( PTT^ ) and yields 



p{x ,x ) 



2 7r2 



(0), 



(1), 



(2), 



+ . . . 



(2.24) 



(2.25) 



Let us for simplicity concentrate on this component of the energy-momentum tensor: all others can 
be also easily constructed using Eq. (|2.12| ). 

Following the original suggestion of [51] (see also [12]) let us first consider delta-function shock 



waves 



ti(x ) 



fii6{x-), t2{x+) = i22S{x+). (2.26) 

As was argued in [1], the delta-function shock waves give a solution of Einstein equations having 
correct qualitative features of the solution for any colliding shock waves with non-negative ti and 
t2- However, plugging the delta- functions from ( p.26|) into Eq. ( |2.16|) and then into Eq. (|2.19|) we 
immediately encounter a problem: we obtain products of delta-functions and theta functions, like 
S{x^) 6{x^). To properly handle those terms let us regulate the delta-functions by spearing them 
along the light cone directions: 



ai 



[X 



-)e{ai-x-), t2{x+) 



^9{x+)e{a2 
a2 



X 



(2.27) 
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To be more specific let us consider in the boundary theory a coUision of two ultrarelativistic nuclei 
with large light-cone momenta per nucleon , P2 , and atomic numbers Ai and A2. In order to avoid 
suppression in each graviton exchange coming from the Newton's constant (see e.g. Eq. ( |2.3| )) 

let us assume that each nucleon in the nucleus has A^^ nucleons in it. This factor of A^^ (^1 ) 

and (T2++) cancels the factor of in Eqs. ( |2.3| ) and ( p.4| ). In the end one obtains, similar 

to [1,73] 

/il~p+A2A}/^ /i2~P2-A^A^/', (2.28) 

while the Lorentz-contracted widths of the nuclei are 

^1/3 ^1/3 

ai i^, a2 (2.29) 

Pi P2 

The scales Ai and A2 are the typical transverse momentum scales describing the two nuclei [1], 
similar to the saturation scales. 

Using Eq. ( ^1271) along with Eqs. (^[25]), ( pJfSD , {^J^, ( p:22|) , and ( p:23a| ), we obtain 



Ar2 

p(x+,x") = -^8nii22X^ X' e{x^)e{x') \l - 12/iix+ {x'f - 12 /i2 {x^fx' + ...]. (2.30) 

2 TT^ 



In arriving at Eq. ( p.30| ) we neglected terms suppressed by powers of ai/x and a2/x~^. Thus 



Eq. ( p.30| ) is only valid when 



— < 1, ^ < 1. (2.31) 

However this is not the only constraint on applicability of Eq. ( p. 301) : requiring that o(ai/x~, 02/0;"'") 
corrections to the NLO terms are much smaller than LO terms, employing Eqs. ( p.28|) and (|2.29| ), 



and assuming for simplicity that ~ P2 > ^1 ~ ^2 = A, Ai ^ A2 = A, we obtain another 
restriction 



AA^/^r<l (2.32) 



with the proper time r = v2x+x^. Hence Eq. (|2.30|) is valid at relatively early proper times 
and acquires order-1 corrections at later times. Indeed Eq. ( p.30| ) provides an exact solution in 
the formal limit of ai, 02 — which reduces ti and t2 back to the delta-function expressions given 
in Eq. ( p.26| ). However, Eqs. ( p.28| ) and ( |2.29| ) demonstrate that if we keep track of the physical 
origin of the delta-functions, the infinitely-thin nucleus limit gets more involved. Of course one can 
always postulate the nuclei to be very thin in the longitudinal direction while keeping their atomic 
numbers fixed, thus making the formal ai, 02 — limit possible: such limit is not attainable in real 
life, but it is a mathematically well-defined procedure. 

Eq. ( p. 30 ) agrees with the appropriate result obtained in [12] for delta-function shock waves. 
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2.2.3 Delta-Prime Shock Waves at NLO 

In [1] it was argued that delta-function shock waves considered in Sect. p.2.2| come to a complete stop 
shortly after the collision, possibly leading to a formation of a black hole. For the boundary theory 
this implied that the colliding nuclei stop after the collision and thermalize leading to Landau-like 
hydrodynamics [57]. This scenario would lead to strong baryon stopping in the collisions, which 
is not what is observed by the experiments at RHIC. Combined with the many successes of small- 
coupling based approaches in describing RHIC data sensitive to early-time dynamics (for a review 
see [50]), this led us to conclude that one can not adequately describe entire heavy ion collision 
within a strong coupling framework. Thus collisions of delta-function shock waves in AdSs are not 
relevant for the heavy ion collisions, in which it is very likely that the initial stages of the collisions 
are weakly-coupled. In [1] to try to mimic these weak coupling effects we suggested using unphysical 
delta-prime shock waves 

t^{x-) = AlAy^6'{x'), t2{x+) = Al Al^^6'{x+). (2.33) 

The shock waves in Eq. ( p.33| ) are fundamentally different from those in Sect. 2.2.2| as the integrals 
of these shock wave profiles over all x~'s and/or x"'"'s give zero. The shock waves ( |2.33| ) have 
unphysical energy- density on the light cone. However, in the LO calculations carried out in [1] it 
was shown that the behavior of the produced matter in the forward light cone of a collision of two 
shock waves ( p.33|) gives a well-behaved physical distribution of matter. This should be contrasted 
with the physical shock waves in Sect. |2.2.2| , for which, due to nuclear stopping, the remnants of 
the colliding nuclei would deviate from their initial light cone trajectories and drift into the forward 
light cone. 

To use the shock waves of Eq. ( |2.33|) for calculating the NLO contribution to the transverse 
pressure p we have to regulate them. We do that by rewriting (|2.33|) as [35, 38] 



,1/3 



ti{x' 



t2(x+) 



A ' 

Al ± 6'i 



X 



-^i ) 



i=l 



(2.34) 



i=l 



Each delta-prime in Eq. ( |2.34D corresponds to a thin slice of a shock wave (a "nucleon" ) localized 
around the longitudinal coordinate xf. The coordinates x~ are localized to the interval [0, ai] of 
the x~ axis, while the coordinate xf are localized to the interval [0, 02] of the x~^ axis. 

Employing Eq. in Eqs. ( ^1251) , {^A§, (^) , and j ^^ ), and assuming that 

Ai,A2 ^ 1, yields the transverse pressure 



p{x^ 



.x-) = f,8AlArAlAr 



[X 



[X 



1 -40 



-36 



AlA\/' 

Al A^' pt x+ {x^ + Al P2 {^^f^ 



1/3 



(2.35) 
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Eq. ( p.35| ) is derived in Appendix ^ Just like with Eq. ( p.30 ), in arriving at Eq. ( |2.35|) we have 
neglected terms suppressed by additional powers of energy, i.e., we assumed the condition ( p.31|) to 
be valid. At the same time we did not have to assume that the bound (|2.32|) applies. 



From Eq. (|2.35|) we see that NLO corrections in the transverse pressure are of two types: 
they can be rapidity/energy-independent, like the second term in the square brackets, which is 
proportional to x~ ~ r^. They can also be rapidity/energy-dependent, like the last term in the 
square brackets in Eq. ( p.35| ), which is proportional to, say, ~ r'^e^, where we defined 

the space-time rapidity t] = (1/2) ln(x+/x^). That term also includes explicit powers of the large 
momentum components Pi and p^, i.e., it is explicitly energy-dependent. Indeed if x~ ^ 1 or 
P2 x~^ ^ 1 the last term in the square brackets of Eq. ( |2.35| ) dominates over the second term in the 
brackets. 



2.3 NNLO Results 

Evaluation of the NNLO terms goes along the same lines as the NLO calculation. One plugs the 
expansion of Eq. ( p.l4| ) into Einstein equations (|2.5|) and expands the resulting equations up to the 
quadric order in t's. In particular one obtains the following equations for G^"^^ and H^'^^ 



Cf) + 5 Hi') _ , Hi') + 2 z i/f ) _ +Ajz' + Ag + An z'' + A,, z'' + A15 z'' = 



(2.36a) 



[zz) Gf^ + 2 Hi'^ - z Gf} -2z Hf} + 48 A-j / + 80 z^ + 120 2" + 168 ^13 2" 

+22AAi^z^^ = 

(2.36b) 



with A's and A's being some known functions of ti and ^2- Eliminating G^'^ from Eqs. ( |2.36a|) and 
(136^) yields 



-3 Hi'^ + 3z Hi'] - z' Hi'], + 2z' H^^}^_ , + 12^rz^ + 16 ^9 / + 20 z^^ + 24 ^13 z^^ 

+28^i5;z^^ = (2.37) 



with 



tj tl tl 





tl tl 




^2 ^2 t2 t2 t2 t2 

Figure 5: Some of the graviton diagrams corresponding to the NNLO solution of Einstein equations found 



in Section 2.3 



\I/7 
^9 



13 



15 



(2.38a) 



a; ; 1 



^11 = — 



1 

60 



128 h^ti t2 + 34 {hQx+x~Y ~ 13 /?.ox+z+ ^ox-x- 
+ |36t2-f^{t-a;- ~ 10^0x+ ^Ox+x-x" - 6 A2 /^o x+ x" x" x" + (^^^ <^ o:"; 1 ^ 2)| 
768 tl t2 /iox+x- ~ 16 hQx+ x-x- hox+x+x- + |l36 ti t'2 /^ox- 
13/iox+x-x-x- /iox+x+ + (a;^ ^ a;~; 1 ^ 2) 



576 



+320^2^.^^^ 



2 x~ x^ 



368t^t2 ~ ^Ox+x-x-x- ^Ox+x+x+x- + 1^70^2 -f^3^-^- + 19 tl ^Ox+x" 



504 

+ (x'*' ^ a;~; 1 2) 



2)_ 
(2.38b) 

(2.38c) 
(2.38d) 
(2.38e) 



The prime in t'i{x ) and in t^ix^) indicates derivatives with respect to the only argument of the 
functions. 

To find a causal solution of Eq. ( f^.37D one expands i/*^^-* into a series in z^, matches the co- 
efficients of the powers of and requires the series to terminate at some finite order to find the 
coefficients. The answer then reads 



X \ x ,z) 



H^^\x^,x-) + ifr(a;+,x-) z« + E)^\x'^ .x') z« + H'^\x-^ ,x~) z 



r(2) 



(2), 



r(2)/ 



^Ef (x+, X-) z^^ + Hf\x+, X-) z^* (2.39) 



(2), 
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with 



H'^' = ^Hf^ (2.40a) 

^ _ 16 _ (16)(30) _ (16)(30)(48) _ (16) (30) (48) (70) ^ 

(2.40b) 

H^'^- ^ 1- 30 (30) (48) (30) (48) (70) 

1 48 (48) (70) 

^'^^-Wa:)^--(^^- ^^-^^^^ 

i/f (2.40f) 



Using Eqs. ( |2.39| ) with ( p.40[ ) in the remaining Einstein equations allows one to find the other 



components of the metric at the same order: G^'^\ F^'^\ and F^'^\ The essential classes of diagrams 
resummed at NNLO are shown in Fig. |^. They involve either three rescatterings in one nucleus and 
one rescattering in the other nucleus or two rescatterings in each of the nuclei. 

3. Asymmetric Collisions of Shock Waves in AdSs 
3.1 Derivation of the Equations 

We now want to find the solution of the proton-nucleus scattering problem at strong coupling. In 
other words we want to resum all-order graviton exchanges with one shock wave while keeping only 
terms with a single graviton exchange with the second nucleus. That is, we want to resum all powers 
of, say, t2, while keeping only the leading power of ti. An example of a typical diagram which is 
resummed this way is shown in Fig. |^. 

To resum the diagrams of the type shown in Fig. ^ let us first construct the corresponding 
Einstein equations describing this classical graviton field. We start by writing the metric, which is 
just the same as given in Eq. ( |2.7|) , but without capitalizing the unknown functions, to distinguish 
from the case of the full nucleus-nucleus collisions: 

ds'^ = ~ \2 + q(x~^ , x~ , z)] dx~^ dx~ + \u(x') + f(x~^ , x~ , z)] dx~'^ 

+ 



- [2 + g{x^,x~,z)] dx^ dx~ + [ti{x') z^ + f{x^ 
t2(a;+)/ + /(x+,x",z) dx^^ + [l + h{x^,x~,z)] dx']_ + dz"^^. (3.1) 



We now want to plug the metric ( |3.1| ) into the Einstein equations ( |2.5| ) and linearize it in ti. In 



doing so we have to remember that, as /, /, g, and h should have only one factor of ti in them. 
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tl 




Figure 6: A diagram contributing to the metric of an asymmetric collision of two shock waves as considered 
in Section ^. 

one has /, /, /i ~ ti. Thus one has to linearize Einstein equations in ti and in /, /, g and h. The 
relevant equations are 

-4:zH2f -8zUit2- zU2f, +g, + 5K-zh,, + zH2K-^- +2zK+^- =0 (3.2a) 
(zz) 8zH2f + 32zU,t2 + 3zU2f, + g, + 2K + zH2f,,-zg,,-2zh,, = (3.2b) 

{-z) z'^t[t2 +z^t2f^- -^gx-z-h^-z-^fx+z = 0, (3.2c) 

where we suppressed the arguments of all functions and, as usual, the subscripts z, and x~ 
indicate partial derivatives with respect to these variables. Again, the prime in t[{x~) (and in 
^2(3;''') below) indicates a derivative with respect to the only argument of the function. Other 
components of Einstein equations are not needed, as Eqs. ( ^.2] ) contain enough information to find 
/, g and h. In fact we will need to know only h: as was shown in Eqs. ( p^.l2[ ) we can reconstruct 
the whole energy-momentum tensor of the produced matter from it. 

Solving Eq. ( 3.2a ) for gz and using the result to eliminate g from Eq. (|3.2c| ) yields 



1 X d X 

-/i. [hzz + zH2{4zHi + fz-zK-,-) -2h,+ ,-] =0. (3.3) 

Eliminating gz from Eq. ( p.2b| ) and solving the resulting equation for fz we get 

fz = -T~rr [-16z'^tit2 -3hz + 3zhzz- hzzz + Az^t2K-x- + t2 K-x-z + 2z^ h^+x-z] ■ 

(3.4) 
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Applying an operator i9-/i9+ to Eq. ( |3.3| ) and substituting from Eq. ( |3.4|) into it we obtain the 
following equation for h 



-3h^ + 3zh^;,- z^h^^^ + 2z^ K+ x- z = lQzUit2 



+ ZU2 



a. 



77 1 2 ^ 1 g 



(3.5) 



Note that the first line of Eq. ( ^.5] ) is identical to the LO equation (4.10) in [1]. Higher order powers 
of t2 come in through the second line of Eq. (|3.5| ) . 

Eq. ( p3| ) is the equation we need to solve. We shghtly simplify it by writing it as 



Z 



16zUit2 



a; 



--zH2d^Z^d,h 



(3.6) 



Below we will use Eq. ( p.6|) to evaluate the diagram in Fig. ^ in the eikonal approximation, which 
we will define in the next Subsection. 



3.2 Green Function and the Eikonal Approximation 



To construct the solution of Eq. (|3.6| ) we will need to construct the retarded Green function 
of the operator on its left hand side. As inverting z'^ dz is trivial, we will need the function 

G{x^, x~ , z; x'~^, x'~ , z') such that 



-dz-dl + 2 d+ a_ 

z 



G{x+,x-, z; x'+, x'-, z') = 6{x+ - x'+) 6{x- - x'-) 6{z - z'). (3.7) 



This is a bulk-to-bulk scalar field propagator, which has previously been found in [74]. For com- 
pleteness of the presentation let us briefly outline the construction of G[x'^ 

Fourier-transforming Eq. ( |3.7| ) into light-cone momentum space (i.e., going from x~^ and x~ 
coordinates to their conjugates k~^, k~ but keeping the coordinate z) and dropping the delta-function 
on the right one can see that the solution of the resulting equation is simply J2{zV2k^ k^). Using 
these Bessel function and going back to the x^ coordinates instead of x^ we write for the 
retarded Green function 

00 00 ^ , , 

r/™0 ^iO\ . /™2 I L2l 

,3_^/3 



G{x ,x ,z;x ,x ,z ) = / dm / dk 



271 J J yVn^ + p 

1 



gi k (x-'—x''') 



xmz'^ J2{mz) ^ J2{mz'). (3.8) 
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The integral over the momentum variable k can be performed yielding 



G{x'^,x ,z;x'^,x' , z') = -9(x^ — x'^)9{x — x' ) — j dm 



X m Jo \/2 (x+ — {x~ — x'~)j J2{m z) J2{m z'). (3.9) 
Eq. ( p.9|) can be further simplified by integration over m, which gives 

x", 2; x'-, 2') = — - - x'-) eis) 6(2 - s) A, ^ ^ ~ (3.10) 

2 TT z'^ (2 - s) 



with 



2 (x+ — x^) [x — X ) — [z — z 



f\2 



(3.11) 



However Eq. (|3.9|) is really all we need for the calculations to follow 

Eqs. (|3.9|) or (|3.1CI|) give us the propagator of the 
gravitons in the s-channel of Fig. ^. These expressions 
allow us to construct the eikonal approximation for the 
graviton production in asymmetric shock wave colli- 
sions. The space-time structure of graviton production 
in such collisions in shown in Fig. |^ It illustrates the 
diagram in Fig. |^: first the graviton (the wavy line) is 
produced in a collision of the proton shock wave and 
some elements of the nucleus shock wave (a nucleon in 
the nucleus). This generates the LO factor of ti ^2- Sub- 
sequently the graviton rescatters in the nucleus shock 
wave with each rescattering bringing in a factor of ^2- 
After the graviton leaves the shock wave it simply prop- 
agates freely. Indeed the transverse dimensions 
and the 5th dimension in AdSs are implied but not 
shown in Fig. |^. 

Most importantly, the propagation of the gravi- 
ton between two successive rescatterings in the nucleus 
shock wave happens over a very short interval in the 
light-cone "plus" direction. Namely the intervals Ax^'s between the rescatterings in Fig. |^ are 
Lorentz-contracted and are all of the order Ax^ ~ 1/p^. As (along with the comparable scale 
pf) is the largest momentum scale in the problem we conclude that Ax^'s are the shortest dis- 
tance scales in the problem, i.e., they are very small compared to any other distance scale. This 




proton 



nucleus 



Figure 7: The space-time structure of the 
graviton emission in a proton-nucleus colli- 
sion. The graviton is denoted by the wavy 
line. After being produced the graviton 
rescatters in the nucleus and the propagates 
freely in the forward light-cone. 
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is illustrated in Fig. ^ which depicts the propagation of the graviton through the highly Lorentz- 
contracted nucleus. Therefore we can approximate the full s-channel graviton propagator by its 
short-x^-interval version. We will call such approximation an eikonal approximation in analogy 
with the terminology used in high energy scattering in four dimensions. 

Putting ~ x'~^ in Eq. ( p.9| ) we can put Jo(0) = 1 which yields the Green function in the 
eikonal approximation 

oo 

1 _ _ z'^ f 

Geik{x~^,x , z; x'"*" ~ x"*", x' , z') = - 6{x~^ — x'~^) 6{x — x' ) — / dmm J2{m z) J2{m z') 

2 z' J 



= ^ 9{x+ - x'+) 9{x~ - x'') 6{z - z). (3.12) 

For inhomogeneous equations like ( p.6|) , or like the following equation 

3 

-h^ - hzz + 2h^+^- = R{x^,x',z), (3.13) 

in the solution, the Green function acts on some function R{x~^, x~ , z) on the right hand side, such 
that 

oo oo oo 

h{x^,x~,z) = I dx'^ I dx'^ I dz' G{x~^ , x~ , z; x'~^ , x'^ , z') R{x'~^ , x'~ , z'). (3-14) 







Using the eikonal Green function ( p.l2| ) in Eq. ( |3.14| ) yields 



heik{x^,x ,z) = ^ R{x^,x ,z) (3.15) 



with the inverse derivatives defined in Eq. ( |2.13| ) 



Going from Eq. ( p.l3|) to Eq. (|3.15|) clarifies the procedure for the eikonal approximation: 
simply neglecting all 2;-derivatives on the left hand side of Eq. ( p.l3| ) compared to d+ we obtain 
Eq. ( p. 151 ). Indeed (9+ ~ l/Ax"*" ~ if the propagator in question spans a short interval Ax"*". 



Hence the main rule of the eikonal approximation is that 9+ is much larger than any other derivative 
in the problem. The short interval scalar field bulk-to-bulk propagator is 

oo oo oo 

dx'^ J dx' j dz Geik{,x^ ,z]x^ ,x'' ,z') [. . .] = q [•••]■ (3-16) 

— oo — oo 

One has to keep in mind that the eikonal approximation should be applied to short-lived 
propagators only. That is we can not just take Eq. ( p.6|) and drop all terms not containing 
As can be seen from Figs. |^ and |^, the graviton propagator after the interaction with the nucleus 



- 19 - 



is not limited to any short interval in any direction. That is, we have to use the full propagator 
( p.9|) for that line. Note that we are not calculating the graviton production amplitude: we are 
calculating the graviton field. Hence in the diagram in Fig. ^ the outgoing graviton propagator is 
off-mass shell, and is not on mass shell, as it would have been for the production amplitude. (See 
e.g. [36] for an example of constructing Feynman diagrams corresponding to classical fields.) 

The graviton propagator in Eq. ( |3.16D does not take into account rescatterings and only de- 
scribes free propagation for a graviton over a short time interval. The eikonal approximation should 
also be applied to the multi-graviton vertices in Fig. To facilitate the application of the eikonal 
approximation let us recast Eq. ( p.6|) in a slightly different form. Defining 



h = z d^h 



Q_ 

h = 16zUit2 + Z^t2^ 



we rewrite Eq. ( p.6|) as 
Di + 2d+ d_ 
where we have defined differential operators 

Di = z^d, 

and 

D2 = z^d.. 

Defining the truncated amplitude 



Do -2d, d. 



z^ z^ ^ 



h--zHod'ih 



(3.17) 



(3.18) 



h 



(3.19) 



(3.20) 



(3.21) 



allows us to write Eq. (|3.18|) as 



9_ 

h = 16zUit2 + zU2^ 

04- 



The solution of Eq. (g) is 



D2-2d+d^) (Di + 29+9_) '-^zU2d^_ [Di + 2d+d^^ 



h. 
(3.22) 



h 



A D2 
X 16z'^tit2. 



Di 
2d.d. 



Di 

2d.d^ 



(3.23) 
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As ^2 and 1/9+ do not commute, here and throughout the paper we have 

2 



^2 IT 

0+ 



that is, each 1/9+ operator acts on everything to its right. 
A simple algebra gives 

-1 



h 



with 



8,-81 + 28+ 9_ 



(3.24) 



1 fh 



(3.25) 



..]{z)= dz'[..]{z') 



(3.26) 



Therefore and h are related to each other with the help of the Green function ( |3.9| ). Therefore 
h is really the part of the amplitude in Fig. |^ without the last s-channel gluon propagator, i.e., 
h is the truncated amplitude. As all s-channel graviton propagators in the truncated amplitude 
h are short-lived, we can apply the eikonal approximation to the equation ( |3.23|) for h. In fact 
Eq. (|3.23| ) is already cast in the form designed to simplify the expansion in inverse powers of 9+. 
The eikonal h we obtain this way can be used in Eq. (p.25|) with the full Green function ( p.9|) to 
recover h{x^ , z). As Eq. ( p.23|) appears to be prohibitively comphcated to evaluate analytically, 
the eikonal approximation appears to be the only way to proceed. In fact, as we will shortly 
see, for the delta-function shock waves it yields the exact solution for the metric generated in the 
asymmetric (proton-nucleus) collision of two shock waves. 

3.3 Delta-Function Shock Waves 

3.3.1 Energy-Momentum Tensor of the Produced Medium 

Let us again consider a collision of two physical delta-function shock waves with ti and ^2 given by 



Eq. ( p.26|) . We will also keep the smeared shock waves in Eq. ( p.27|) in mind. 

First let us apply the eikonal approximation to Eq. (|3.23| ) without substituting the explicit 
expressions for ti and ^2 from Eq. ( p.26| ). As we argued above, in the eikonal approximation the 
derivative 8+ is the largest momentum scale in the problem. Hence in Eq. (|3.23|) we have 



<1, 



Do 



2 8+8 

After neglecting those terms Eq. (|3.23|) yields 



2 8+8^ 



< 1. 



(3.27) 



eik 



1 



2 8+ 



16zUit2. 



(3.28) 
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(3.29) 



To evaluate Eq. ( |3.28|) we expand it in a series 
Using Eq. in Eq. (^) along with Eq. (^) we write 

oo oo oo _ 

heik = J dx'^ J dx'~ J dz' G{x^ , , z; x'^ , x'~ , z') — (^—j^^ 

— oo — oo 

x~ oo ^ oo 

dx'^ I dx" I dz' -— I dmm Jq (^m a/2 (x+ — x'+) (x^ — x'^)^ J2{m z) J2{m z') 

{z')^''+Hi{x'')t2{x'+). (3.30) 



n=0 



n + 1 / 1^ 
^ 2n + 3 V 2 



9' 



d' 



Here = d/dx'^. The expression ( p.30|) is still rather difficult to evaluate. However, as we are 
primarily interested in the dynamics of the gauge theory, we only need the z'^ term in this expression 
to obtain the transverse pressure of the produced medium using Eqs. (|2.12|) and ( |2.24|) . As the 
series expansion of the Bessel functions converges everywhere, we expand 



J2{m z) = -m z + o{z ) 
8 

in Eq. (|3.30| ) and integrate over z' and m obtaining the transverse pressure 



(3.31) 



P 



Ar2 

c 

2 7r2 



8 ^(n + l)2(-2)" / / dx'-(x+-x'+)i+2"(x--x'-) 

n=0 



-oo — oo 



X {d'!^t,{x'-)\ [t2{x'^)—] t2{x'+). 



(3.32) 



Using integration by parts in the integral over x' in Eq. ( p.32| ) and remembering that ti is a 
localized function of x~ yields the final expression for the transverse pressure 



n: 



2n 



%8 5^(-2)"(r^ + l) 



{2n + l)\ 



n=0 



dx'- (x" -x'-)^+"ti(x'~) 



X j dx'+{x+ -x'+f+^^ [t2{x'+)—] t2{x'+). 



(3.33) 
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Eq. ( 3.33 ) is one of the main results of this Section. It is the simplest expression for p we could find 
without using an explicit form for the functions ti and ^2- 

As with the NLO calculations of Sect. |2.2.2| , substituting ti and ^2 from Eq. ( |2.26| ) into Eq. ( p.33|) 
would generate terms like 5{x'^) 9{x^), evaluation of which is ambiguous. To avoid this ambiguity 
we use the smeared ti and t2 from Eq. ( 2.27|) . For ^ ai and ^ 02 we have 



dx 



'- 'x- 



X 



/il [x 



Q[x- 



(3.34) 



and 



t2(x 



1 



(n + l)\ 



(x+)^+2"^(a;- 



(3.35) 



Using Eqs. ( p.34|) and (|3.35|) in Eq. (|3.33|) and summing the series over n yields 



P 



8fiifi2X^x-9{x+)9{x-) 
2vr2 [1 + 8/X2(x+)2x- 



l3/2 



(3.36) 



This is another main result of this Section: Eq. ( |3.36| ) gives us the transverse pressure of the 
medium produced in the collision of a proton and a nucleus at strong coupling. It is valid at 
x~ ^ ai and x"^ ^ 02: these conditions are automatically satisfied if the sources are exact delta- 
functions of Eq. (|2.26| ). Hence for the delta- function sources ( p.26|) Eq. ( ^.361) provides us with the 
exact solution of the problem! As a cross-check one can see that expanding Eq. (|3.36| ) in a series in 
/i2 to NLO yields Eq. ( CT) - 

Eq. ( p.36| ) allows us to explicitly specify the limits of our approximation. Namely, we resum all 
powers of rescattering in the nucleus, which, for delta-function shock waves translate into powers 



of H2 (x' 



At the same time we neglect higher rescatterings in the proton, which, by analogy, 



would bring in powers of Hi{x' 



x^ 



Hence the apphcability region of Eq. (|3.36|) is defined by 



/il (x )^X"^ ^ 1, IJ,2{x^Yx ~ 



1. 



(3.37) 



(Indeed for small ^2 {x^Y ^ Eq. (|3.36|) applies too.) For non-delta function shock waves like those 



given in Eq. (|2.27|) one also has to keep the limit ( p.32|) in mind while studying the applicability 
region of Eq. (|3.36|) . 

Using Eq. ( p.36| ) along with Eq. (|2.12|) we can find all other non-zero components of the energy- 
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momentum tensor of the produced medium: 

(T^^)^-4^^-^-(-^)^^(-")t'^ (3-38a) 
27r2 [l + 8/i2(s+)2x-]'/' 



2 
3 



X 



3 v^l + 8 At2 + 4 /is (^9 + 16 /i2 - 6 v^l + 8 /ia x") 



[1 + 8/^2 (x+)2x- 



-|3/2 



(3.38b) 



A/'2 8/ii/i2X+x-^(x+)e(x-) 



(r"-) = s : ' ^ (3.38c) 



2 [1 + 8 /i2 (x 



8/ii/i2X+x-6'(x+)6'(x") 

2 71-2 ll + SfX2{x+yX- 



T') = S . , '1,2 ■ (3-38d) 



Provided the complexity of the problem at hand, the resulting formulas (|3.38| ) for the energy- 
momentum tensor are remarkably simple! 

Now we can ask a question: what kind of medium is produced in these strongly coupled proton- 
nucleus collisions? Is it described by ideal hydrodynamics, just like Bjorken hydrodynamics was 
obtained in [51]? In our case the produced matter distribution is obviously rapidity-dependent, 
so it is slightly more tricky to check whether Eqs. ( |3.38| ) constitute an ideal hydrodynamics, i.e., 
whether it can be written as 

T^^ = [t^V)u^u'' -vri^'' (3.39) 

with the positive energy density e and pressure p. ri^^ is the metric of the four- dimensional 
Minkowski space-time and is the fluid 4-velocity. 

For the particular case at hand it is easy to see that the energy-momentum tensor in Eq. ( p.38| ) 



can not be cast in the ideal hydrodynamics form of ( ^.39| ). In the case of ideal hydrodynamics one 
has 

T++ = (e + p) {u^f > 0. (3.40) 

At the same time (T"*""^) in Eq. ( |3.38a| ) is negative deGnite. Therefore the ideal hydrodynamic de- 
scription is not achieved in the proton-nucleus collisions. We believe this result is due to limitations 
of this proton-nucleus approximation. Any strongly coupled medium at asymptotically late times 
is almost certainly bound to thermalize. Our conclusion is then that thermalization/isotropization 
of the medium does not happen in the space-time region defined by the bounds in Eq. (|3.371 ). What 
we found in Eq. (|3.38|) is a medium at some intermediate stage, presumably on its way to ther- 
malization at a later time. It is likely that one needs to solve the full nucleus-nucleus scattering 
problem to all orders, as shown in Fig. 0, to obtain a medium described by ideal hydrodynamics. 
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3.3.2 Proton Stopping 

In [1] it was argued tliat tlie pliysical sliock waves given by Eq. ( |2.26| ) or by Eq. ( p^.27D come to a 
complete stop shortly after the collision. The conclusion was based on the LO calculation, which 
for the shock waves ( p.271 ) gave the following ++ component of the energy momentum tensor of a 
nucleus (or its remnants) moving in the light-cone plus direction after the collision 



(T 



X 3> a, X 



a/2)) 



2 7r2 a 



[l-2/ix+^a] 



(3.41) 



In arriving at Eq. ( |3.41| ) in [1] we for simplicity put /ii = /X2 = and oi = 02 = a. Eq. ( p.41|) allowed 
us to conclude that as the light-cone time ~ a the -|--|- component of the energy momentum 

tensor of the shock wave would become zero, meaning that the shock wave stops propagating along 
the light cone. Indeed, as we saw above (see Eq. ( |3.37] )), at the same time as the stopping happens, 
i.e., when ^x^"^ a ~ 1, higher order graviton exchanges would become important. With the help of 
Eq. ( 3.33| ) we can now explore whether multiple graviton exchanges with the nucleus shock wave 
modify our conclusion about proton stopping reached in [1] at the LO level. 

We start by evaluating Eq. ( p. 331) for x^ 3> a2, but with < x~ < Oi. That way we follow the 
proton shock wave for some time after the collision, which allows us to find the energy-momentum 
tensor of the shock wave itself. As x+ ^ 02 still, Eq. (|3.35|) remains unchanged. We have to 
re-evaluate the left-hand-side of Eq. (|3.34]) for < < ai. This can be readily done yielding 



dx'- (x- - x'-)^+" ti{x'') = ^ (x~)'+", for < X- < ai. 

CL ]^ 7^ I ^ 



(3.42) 



Using Eqs. ( p.35| ) and (|3.42| ) in Eq. ( |3.33|) and resumming the series one obtains the transverse 
pressure inside the proton shock wave, which, with the help of Eq. (|2.12| ), gives the following 
expression for the +-|- component of the energy momentum tensor of the produced matter 



Nl /ii 



- <-l + 



1 



l + 8/i2 (x^ 



for < X < Oi. 



(3.43) 



X 



Eqs. (|2.3| ) and ( 2.27| ) give the energy-momentum tensor of the original incoming shock wave itself 
as 



iV2 



— , for < X < ai. 



(3.44) 



2 7r2 ai 

Adding the energy-momentum tensors of the original shock wave and the produced matter given in 
Eqs. ( p.44| ) and ( |3.43| ) together we obtain the total ++ component of the energy-momentum tensor 
of the proton shock wave 



27r2 ai ^l + 8/i2 {x+f : 



for < X < ai. 



(3.45) 
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Expanding Eq. (|3.45| ) in the powers of /i2 at x = ai/2 would yield Eq. ( |3.41|) , providing an 



independent consistency check. 

Eq. (|3.45|) clearly demonstrates that the ++ component of the energy-momentum tensor of 



the proton shock wave is positive deGnite. Notice that the LO solution ( 3.41 ) for becomes 
negative at large enough x~^. Inclusion of multiple graviton exchanges fixes this problem. T"*""*" in 
Eq. ( |3.45| ) goes to zero smoothly as grows large for any fixed x~ in the < < ai range. Thus 
Eq. ( 3.45 ) explicitly demonstrates that strong-coupling interactions of the proton shock wave with 
the nucleus shock wave would stop the proton shock wave shortly after the collision. For x~ = ai/2 
the stopping happens at x+ ~ 1/^/^2 ai, in agreement with the arguments of [1]. 

3.4 Delta-Prime Shock Waves 
3.4.1 Deltology 

The eikonal approximation used in Sect. |3]^ reduces the exact formula (|3.23| ) to Eq. (|3.28| ). 



Eq. (|3.28| ) resums the powers of t2 with only one factor of 1/(9+ inserted between each pair of 
t2^s. It thus resums terms consisting of sequences like 



see also Eq. (|3.33|) ). This is indeed natural in the eikonal approximation, as (9+ ~ P2 is large and 



we want to have as little powers of 1/9+ ~ 1/^2 possible in each term. Eq. ( |3.28| ) resums the 
absolute minimum number of the powers of 1/9+ . 

An attentive reader might have noticed that the approximation of Eq. ( p.28|) is insufficient for 
the delta-prime shock waves of Eq. ( p.33|) . Indeed performing the calculation in Appendix ^ we 
saw that leading powers of arose not only from the terms of the type shown in Eq. ( p.46| ), like 
we had in Eq. ( [A10| ), but also from terms with two powers of 1/9+ inserted between two t2's, as 
can be seen from Eqs. (|A6| ) and (|A9|). Hence we need to rethink our power counting if we want to 
resum the leading eikonal terms for the delta-prime shock waves. 

Let us start with delta-function shock waves with t2 ~ 6{x^). In this case 

t2 7^ t2 ~ 5(x+) 9{x+) ~ 5(x+). (3.47) 

Here we are being rather sloppy in treating S{x^) 6{x^): of course the whole regularization intro- 
duced in Eq. ( |2.27|) above was designed to obtain the correct values for 9{0) in different situations. 



However, for the purposes of counting powers of P2 the exact value of 6^(0) is not important as long 
as it is a P2 -independent number. Eq. ( |3.47| ) demonstrates that for t2 ~ 6{x~^) one has 
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Therefore t2 ~ o(l) in power counting. 

Higher order corrections may come in through an insertion of one power of l/d"^ between two 
t2^s. One then gets 



h^t2 S{x+)x+9{x+) = 0. (3.49) 



The equahty in Eq. ( p.49|) is only true for delta-function shock waves and demonstrates that 



Eq. ( |3.38| ) is the exact solution for the problem of the collision of two delta-function shock waves. 
For the smeared shock waves of Eq. ( p.27| ) the zero in Eq. ( |3.49| ) would be replaces by 02 S{x~^). As 
a2 ~ ^/P2 this indicates suppression by a power of I/P2 compared to the leading-order terms in 
Eq. ( p.48|) . One can similarly show that insertions of higher powers of 1/(9+ would bring in further 



suppression. Thus our approximation in Sect. |3]^ is justified by this explicit power counting. 



Let us now turn our attention to delta-prime shock waves of Eq. ( |2.33| ). Notice that ti and t2 
in Eq. ( p.33| ) do not explicitly depend on and p^: as we show in Appendix the dependence 
on these momenta (and, hence, on the center-of-mass energy of the collision) comes in through 
singularities like 6{x^ = 0). For t2{x~^) ~ 6'{x^) one has 

t2^t2 ~ 5'(.T+)<5(,T+) ~ ~ P2 5'{x^). (3.50) 

Again we are not keeping track of factors not containing p^. Iterating the procedure we get 

^2^) t2 ~ {P2rS'{x+), (3.51) 

that is, for delta-prime shock waves ^2(1/^^+) ~ P2 ■ 

To understand higher order terms with more powers of 1/9+ consider 

h^t2 ~ 6'{x+)e{x+) ~ {6{x+) e{x+)y - 6\x+) ~ ((5(x+)^(x+))'-j92 5(x+). (3.52) 

The term which was sub leading for the delta-function shock waves (see Eq. (|3.49|) ) gives a leading- 
order factor of P2 for delta-prime shock waves, as we see from the last term in Eq. (|3.52| ). Applying 
higher powers of t2 (1/(9+) to the last term in Eq. (p.52 ) does not make the term less important: 



(t2j^yP2S{x^) ~ (p2-)"+^5(x+). (3.53) 

(In fact the {S{x^) 9{x^))' term in Eq. (|3.52|) also brings in powers of after the operator t2 (l/f^+) 
acts on it at least once.) We thus see that one insertion of ^2(1/^^+) still gives leading terms in the 
case of delta-prime shock waves. Fortunately higher order insertions of t2{^-/d\) start generating 
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subleading terms and can be discarded. We illustrate this by acting with ^2(1/'^+) the last term 
in Eq. (|33^ ): 



No extra powers of P2 is generated and hence such terms are subleading. 

Insertions of a higher number of inverse derivatives are also subleading. For instance 

t2^h 5'{x+) x+e{x+) ~ -(5(a;+)^(a;+), 



(3.54) 



(3.55) 



again producing no powers of p^. 

We conclude that for delta-prime shock waves the eikonal approximation consists of the term 
in Eq. ( |3.28| ) along with all terms with a single insertion of ^2(1/^+) i'^ possible positions. 



3.4.2 Energy-Momentum Tensor of the Produced Medium 

To take into account all leading terms for the delta-prime shock waves we write 



(3.56) 



with /leifc given by Eq. ( |3.28|) and 5h denoting the sum of all terms with all-orders of t2(l/<9+) and 
exactly one insertion of t2{l/d\) as contained in Eq. (|3.23|) . 

Expanding Eq. ( p.23|) to the first order in t2{l/d\) we write 



5h 



. 1 4 9- 
2 ' 



n -1 



^Z''t2^ [Di+D2]+lzU2^zH2^D 



heik- 



(3.57) 



Expanding the first factor on the right hand side of Eq. (|3.57 ) into a series and using the series 
representation for heik from Eq. ( |3.29| ) yields 



5/1 = ^ (m + 1) 



m=0 



2 (9, 



1 



9_ 



1 



-zH2^ [Di + D2) + -zH2 — zH2^ 



J2{n + 1) (--zH2-^] 16zUit2. (3.58) 



Using the definitions of Di and D2 from Eqs. (|3.19| ) and ( |3.20| ) we obtain 

^4n+7 _ _g^n+l){2n + 1) 2^"+^ 

and 

+ D2) = -4 (n + 1) (2 n + 3) /"+^ 



(3.59) 



(3.60) 
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which allow us to rewrite Eq. ( p.58| ) as 

oo 

6h = 16 (n + l)^(m + 1) 



n,m=0 



n+m 



-2(2n + 3) z 



4 (n+m)+9 



t2 



dl 



oi 



{^2^ h. (3.61) 



The truncated amplitude contribution (|3.61| ) leads to the contribution to the amplitude through 
Eq. ( p. 251) . Using the Green function ( p.9[ ) we write 



oo oo 



6h 







2^2 y dx'^ J dx'' j dz' j dq q Jq ^2 {x+ - x'+) {x- - x'-)^ J2{q z) J2{q z') 

oo / 1 \ / 1 \ 

{n + lf{m + l) -- ^2(x'+)^ 



2 (2n + 3) ^^,^^4(„+„)+7^^^__^,+ -^ 2n+l ^^/^|4(„+m)+ll ^, 



n + m + 2 



df n + m + 3 

X [9'"+™ti(x'-)] (t2{x'^)l-j t2(x'+). 

(3.62) 



As in Sect. p.3| we are interested in the contribution of the metric element, now the term Sh, to 
the transverse pressure. As further evaluation of Eq. (|3.62| ) to all orders in z appears to be rather 
involved, we expand it to the order z'^ using Eq. ( |3.31| ), integrate over z' and q, and eliminate the 
d'_ derivatives by successive integrations by parts. This yields the following contribution to the 
transverse pressure 



6p 



Tvr2 r r ^ 

/ dx'+ / dx'- J2 (^ + l)'(m+l)(-2)"+'"(x- -x'-)"+'"+2^i(x'-) 



— oo — oo 



n,m=0 



X (x+ - a;'+)2{n+m)+2 

+ 4(2n+ 1 



2(2n + 3) 



(2 (ri + m) + 3)! 
{n + m + 2)! 



(2 (n + m) + 5)! , , , , , 2 , /-x ^ / /+\ 1 

(x+-x+)''(x -X )t2{x^)^ 



(n + m + 3)! 

X {^^^^'^) IrJ (^2(x'+)^)%2(x'+). 



(3.63) 



For further evaluation of Eq. ( p. 631) we will explicitly substitute the shock wave profiles from 
Eq. (|2.33| ). As one can see from the calculations in Appendix ^ the regularization in Eq. ( p.34|) is 
not needed for the leading-p^ terms. 
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First let us find the contribution to the transverse pressure coming from the piece in Eq. ( 3.33 ), 
which we will refer to as Peik- As can be easily shown for ^2 from Eq. (|2.33| ) (see Appendix 0) 



t2(x+) 



^4 



t2{X j - -, — TT 

[n + 1)! 



(3.64) 



In Eq. ( |3.64| ) and henceforth for simplicity we absorb factors of A^^ and A^^'^ into A^ and A^. Using 
Eqs. (^ and (^) in Eq. yields 



Peik 



iV2 

c 

2^ 



8Al kle{x+)e{x- 



1 -44p2 A^(x+)2x- + 64 [p^ Kl{x+f x~y 



[1 + 8P2-A2 



[X^ 



' X' 



7/2 



(3.65) 



To evaluate Eq. ( 3.63|) one needs another relation (see Appendix for its derivation), valid 
only for ^2 from Eq. (|2.33|) at the leading order in 



t2(x+) 



1 



1 

a: 



-1)"^+^ (A2)"+-+2 (p- 

{n + 1)! (m + 1)! 



\n+m+l 



5(x+). (3.66) 



Using Eqs. ( ^1331) and in Eq. we get 

iV2 



5p 



27r2 



8A2 A2^(x+)^(a;-) 



X <^ 1 -36j9^A^(x+)2x" 



1 -44p^A2(x+)2x- + 64 (p2 A|(a;+)2a;-)^ 



[l + 8p2 Ai {x+Yx-]"^ 



(3.67) 



The net transverse pressure is obtained by adding Eqs. ( p.65| ) and ( |3.67|) 



P = Peik + Sp 



2 7r- 



8A^A2^(x+)0(x') [1-36j9^A^ 



[X 



X 



-]■ 



(3.68) 



This is the first main result of this Section. Importantly all higher order terms cancel leaving us 
with the simple expression ( |3.68| )! Note that Eq. ( |3.68[ ) is just a sum of the LO and NLO corrections 
(resumming leading powers of P2 ); ^nd thus it agrees with Eq. ( |2.35 ). Namely it turns out that the 
NLO transverse pressure from Eq. (|2.35| ) taken at the leading-p^ accuracy gives us the full eikonal 
result for the proton-nucleus scattering problem with delta-prime shock waves. 

Eq. ( p.68|) indeed has a limited region of applicability. As it was derived for the proton-nucleus 
approximation, similar to Eq. ( p.37|) we must have 



p^AKx'Yx^ < 1, p^A^iix^fx 



(3.69) 



to be able to neglect eikonal graviton exchanges with the proton shock wave. We also want to 
neglect the non-eikonal terms shown in Eq. (|2.35D, which requires (see Eq. (|2.32|)) 



AIt^ < 1, Alr^ < 1. 



(3.70) 
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It appears that the region of apphcabihty of Eq. ( |3.68| ) is indeed somewhat hmited and is 
confined to the region of large and small x~ , i.e., the region of space-time in the forward light 
cone bordering the proton shock wave. Still it is perhaps surprising to see that the pressure in 
Eq. ( p. 68 ) can easily become negative at large enough x^. As the pressure is negative we conclude 



that the system has not yet reached the ideal hydrodynamics state, similar to what happened to 
the delta-function shock waves in Sect. |3.3| . 

The presence of negative pressure does not pose any problems by itself: negative pressure is 
known to arise in the early stages of heavy ion collisions when they are described in the Color 
Glass Condensate framework [24,25]. Even in the strongly-coupled theory considered here, the LO 
part of the transverse pressure (|2.35|) , when used in Eq. (|2.12| ), leads to negative pressure in the 
longitudinal direction [1]. In comparison, appearance of a negative energy density would be indeed 
worrisome and would indicate an unphysical situation. However, we can not calculate the energy 
density here, as the matter distribution is indeed rapidity-dependent and is not described by the 
ideal hydrodynamics: it is impossible to find the local rest frame of such medium to meaningfully 
talk about the energy density. 

Therefore negative pressure in Eq. ( p.68| ) may be physical. One could interpret it as follows: 
when we chose the delta-prime shock waves of Eq. ( p^.33D we "forced" the shock waves not to stop and 
to continue along the light cone trajectories. At the same time the produced strongly-interacting 
medium is still trying to pull them back together. As the shock waves are "artificially" pinned down 
to their light cones they do not stop, thus creating a negative pressure in the medium which tries 
to slow them down. 

Using Eqs. (p.68| ), ( |2.24]) and ( |2.12| ) we construct all non-zero components of the energy- 



momentum tensor of the produced medium: 

/T++) = — ^ [-8 A? A^ 6(x~) x+ e(x+) + 96 po Al At (x+f e(x+) e(x-)] , (3.71a) 
2 vr^ '- -' 

{T—) = \-8AjAl5ix+)x-e{x-) + 288p^AlAtx+9{x+){x-f9{x-)], (3.71b) 
2 vr^ 

= 7ri8AlAleix-')e{x-) [l-3Qp^ AUx+Yx'], (3.71c) 

(T'l) = 5^ ^sAlAl e{x+) e{x') \1 - 36 Al {x+f x"! . (3.71d) 

2 TT^ 

This is the second main result of this Section. 

3.5 Validity Range of the Perturbative Expansion 

Before proceeding to the conclusions, let us check the validity range of the perturbative approach for 
solving Einstein equations followed throughout the paper and outlined in Eq. ( p.l4| ). For the sake of 
simplicity, we will explicitly analyze the relative contribution to the _L_L metric coefficient obtained 
at LO (i/W), and NLO (H^^^), for the case of (physical) delta function shock waves of Eq. (|2.26| ). 
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Starting from Eq. ( p.l5 ) for the LO results and Eqs. ( |2.22| - p^23D for the NLO one and dropping 



some trivial factors of order one, the ratio Rnlo/lo between the NLO and LO contribution to the 
metric coefficient H{x~^,x~,z) is given by 

_ gW(a;+,x-,z) 3. 1 + i^M' + {z/ry + (^/r)^ 
^^^^/^^ = if(0)(x+,x-,z) ^ ^ " 1 + {z/ry ■ ^^-^^^ 

For definitiveness we have chosen to concentrate on the higher order corrections due to graviton 
exchanges with the shock wave described by the energy scale /X2- To obtain an estimate for the 
NLO graviton exchanges in the other shock wave one simply has to replace ^2 /^i and 77 — > —r] 
in Eq. (|3.72| ). In arriving at Eq. ( |3.72| ) we have made use of the relation = r e''/V2, with 
1] = (1/2) ln(x+/x~) the space-time rapidity. Note that all the coefficients in the numerator and 
denominator in the ratio in Eq. (|3.72|) are put to be equal to 1 for simplicity of the parametric 
estimate we are performing. 

Similar to Eq. one can build the ratio of the NNLO contribution, H^'^^ given by Eqs. 




Figure 8: Schematic plot of the validity range of our solution in the z-t plane. The z-axis was moved to 
a slight left of r = line for illustrative purposes only: indeed r > 0. The darker shaded area indicated 
the validity region of the perturbative expansion of Sect. The lighter shaded area (together with the 
darker one) depicts the broader validity region of the pA approximation developed in this Section. 
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( p.3S - 2!39|) , to the NLO one again for the delta-function shock waves: 

Kn NLO /NLO = H(.'^){x+,X-,Z) ~ 1" 6 ^NLO/LO- (3.73) 

One can easily show that Eqs. ( |3.72|) and ( p.73| ) are also valid for delta-prime shock waves of 
Eq. (p.33|) . Thus, the condition Rnlo/lo ^ 1 sets the validity range of the whole perturbative 
expansion. Two different situations can be considered: [i) z <^ t: in this case the perturbative 
expansion is justified as long as r < [l/(/i2 e'')]^/^. {ii) z ^ r: in this case Rnlo/lo ^ 1 if 
^ ~ [''"/ (/^2 e^)Y^^- The validity region of our perturbative expansion given by the union of regions 
(z) and (ii) is depicted in Fig. |] (see the darker shaded region there). Thus our approximation is 
valid in comparable intervals in r and z near the boundary of the AdS space. 

On the other hand, in the calculation performed in this Section for asymmetric pA collisions, 
we resummed all-order graviton exchanges with the nucleus shock wave. Therefore the applicability 
of this approximation does not require Rnlo/lo ^ 1 for higher order corrections bringing in powers 
of fi2- We do neglect all higher-order graviton exchanges with the proton shock wave that bring in 
higher powers of fii. Thus our pA approximation is valid only if Rnlo/lo ^ 1 for NLO corrections 
with fi2 fJ'i and rj —>■ —rj in Eq. p.72| ). The two regions (i) and {ii) become: (z) z -C r, 
r < [l/{fiie^'^)]^^^, and (ii) z ^ t, z < [t / {^ie~'^)Y^'^ . These new regions (z) and {ii) for pA 
collisions are shown in Fig. ^ by the lightly shaded area. Indeed the validity of the pA approximation 
is much broader than that of the validity of the perturbation series of Sect. |[ 

Importantly, the stopping time estimated in Sect. |3.3.2| , Tstop ~ [Vl/^ae'')]^''^, lies within the 
validity range of our approximation, as one can easily see from Fig. Our approach is valid for a 
comparably broad range of z's at the stopping time, though for the physics of the gauge theory in 
four dimensions only the knowledge of the metric in the small-z region is needed. 



4. Conclusions 

Let us summarize our main results. In Sect. H we constructed the NLO and NNLO terms in the 
perturbative expansion in graviton exchanges for the collision of two shock waves. Our expansion 
generalizes similar expansion constructed previously in [12] from the delta-function-only case con- 
sidered in [12] to the case of shock waves of arbitrary profile. In particular we see that even the 
delta-prime shock waves, which at LO gave rapidity- independent distribution of matter [1], lead to 
rapidity-dependent energy-momentum tensor at NLO (see Eq. ( ^.35| )). 

It is worthwhile noting that the perturbative graviton expansion of Sect. ^ can be built con- 
sistently without introducing a dilaton field. Hence at all orders of the solution the dilaton field is 
zero. As the dilaton field is dual to the (tiF^^) operator in the gauge theory at the boundary we 
conclude that for the collisions of shock waves considered here 

(trF^^) = (4.1) 
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at all times. Thus electric and magnetic modes are always equilibrated in this strongly-coupled 
medium. This result should be contrasted with that of [53], where dilaton field was needed to 
construct singularity-free pre-asymptotics to the Bjorken hydrodynamics metric of [51]. As our 
calculations show, the absence of a dilaton field in the initial shock waves leads to no dilaton field 
throughout the collision. As it is difficult to construct shock waves with non-zero dilaton field it is 
not clear how to construct shock wave collisions with non-zero dilaton field in the forward light cone. 
Therefore the no-dilaton aspect of our result may give one reasons to worry whether dual-Bjorken 
geometry of [51] is obtainable at all in collisions of AdS shock waves. 

In Sect. ^ we have devised an eikonal resummation procedure, which resumed all graviton 
rescattering in one nucleus while keeping only one graviton exchange with another nucleus. The 
results for delta- function shock waves are given in Eq. ( |3.38| ). As is clear from Eq. ( |3.38| ) the 
matter distribution obtained in the proton-nucleus approximation can not be described by ideal 
hydrodynamics, and should be viewed as some intermediate stage of the matter evolution towards 
isotropization. We also showed explicitly in Sect. |^ that strong interactions with the nucleus stop 
the proton completely, as can be seen from Eq. ( p.45| ). 

The eikonal expansion for an asymmetric collision of two delta-prime shock waves terminates 
at the level of two graviton-exchange with the nucleus. The results of the resummation are shown 
in Eq. (|3.71|) . It is important to note that the energy- momentum tensors for delta-function shock 
waves (|3.38| ) and for delta-primes ( 3.71 ) are strongly rapidity-dependent. It is unlikely that a 
matter distribution which is strongly rapidity-dependent at early times would become rapidity- 
independent at late times: such behavior would be acausal, as different rapidity regions become 
causally disconnected from each other as the collision evolves. It is therefore probable that collisions 
of shock waves in AdS will lead to a rapidity-dependent final state at late times: if, due to strong 
interactions, the matter in this late-time state would be described by the ideal hydrodynamics, 
this hydrodynamic description can not be that of rapidity-independent Bjorken hydrodynamics 
[30]. While a full (possibly numerical) study of the nucleus- nucleus collision in AdS would provide 
definitive answer to this question, it is possible that Bjorken geometry gives a good approximation 
to the dynamics of the matter produced in a collision of two identical nuclei only in a narrow interval 
around mid-rapidity. Our results here can serve as a benchmark for further (possibly numerical) 
studies of the collision of two shock waves beyond the asymmetric approximation done here. 

It may be that to obtain Bjorken hydrodynamics in a broader rapidity range one has to abandon 
the idea of colliding shock waves in AdS, and try to simulate the initial stage of the medium by 
matching the AdS metric onto the results for the energy-momentum tensor known from weak- 
coupling CGC methods [24-29]. Such approach was advocated in [14,28,75] and may prove to be 
quite fruitful. One possible shortcoming of the matching method is in the fact that it leaves too 
much freedom in the choice of the early-time AdS metric, leading to a possible loss of uniqueness 
in the description of the subsequent time-evolution of the system. Further research is needed to 
understand which AdS approach is better suited to describe heavy ion collisions. 
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A. Transverse Pressure for Delta-Primes at NLO 



We want to find the NLO contribution to the transverse pressure due to delta-prime sources. Using 
Eq. ( p.25|) we write 
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where we used H^^^ given by Eq. ( |2.23a] ). Our goal is to evaluate the right hand side of Eq. 
for ti and t2 given by Eq. ( |2.34[ ). For simplicity we will evaluate the terms with one power of ti 
only: the remaining terms with only one power of t2 can be obtained by the substitution ti ^ t2- 
Using Eqs. ( |2.19|) along with Eqs. ( p.l6| ) we write 
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for X ^ ai. Similarly we write 
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where in the last step we have used x^ ^ a^. Combining Eqs. (A3) and (A4) we obtain 
6 



(5+5- 



^■j ^ -32 A\ ^J/' (A^ A)!^ x+ X- e{x+) e{x~) + (1 ^ 2). 



(A5) 



- 35 - 



Eqs. (g;T9D, ( pJ6D give 
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The only difference of Eq. (|A6|) with Eq. (^^) t2-dependent part, which is evaluated to give 
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The last term in the last line of Eq. ( [A7| ) is only non-zero when xf = xf , which is only true if 
i = j as all the xf^s are different. However, if i = j that term becomes a delta-function squared. 
Regulating the infinity by the largest momentum scale in the problem we replace 6{x~^ = ^) ~^ P2 
and get 
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for x^ ^ a2- Combining Eqs. ( [A^ ) and (|A8| ) in Eq. ( [A^ ) we get 
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(Indeed the delta function 6{x = 0) should be regulated hj pf.] 
Similar to the above one gets 
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Using Eqs. 
given in Eq. (|2.35|) . 



and ( |A10|) in Eq. (|A1| ) yields the NLO contribution to the pressure 
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B. Iterations of Delta-Primes 

We start by proving Eq. ( p.64| ) for 

t2(x+) = kl5\x+). (Bl) 
Evaluating one iteration of )f:2(l/^+) operator acting on ^2 we get 

t2(x+) t2(x+) = (A^)^ 5\x^) = (A^)^ ' . (B2) 

Similarly 

t2(x+) i-) ' t2(x+) = (A^)^ 6\x^) ^ = (A^)^ ' . (B3) 

It is now straightforward to see what happens at each step of application of the t2(l/<9+) operator 
to write 

which is exactly Eq. (|3.64| ), as desired. 

Now let us prove Eq. ( |3.66|) for t2 from Eq. (|Bl|). First of all use Eq. ( pi]) that we just proved 
to write 
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In each step in Eq. ( p7| ) we neglected a total derivative: it can be shown that those total derivatives 
do not generate leading powers of p^. Eq. ( pTf ) gives 
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where we again neglected the total derivative as it is subleading. In Eq. (|E^ ) we have also regularized 
the extra powers of 6{x^ = 0) by replacing them with p^. 
Combining Eqs. (|B3) and (|E5|) yields 



1 \ 1 
which is exactly Eq. (|3.66| ). 
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